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Abstract 

J> We discuss the prospects for production and detection of fermiophobic bosons (exotic bosons 

decaying only to standard model gauge bosons) at the LHC, and describe simple methods for 
determining spin. We consider two complementary approaches to spin determination: the search 
for decays in the diphoton channel, and the comparison of events with no extra spectator jets to 
<3\ those with one extra jet. We show that these approaches can together allow the fermiophobic 

^ boson's spin to be determined over a wide range of parameter space. We study both even and odd 

t-H parity states. 

> 
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I. INTRODUCTION 



In this paper we will study the phenomenology of massive fermiophobic bosons (denoted 
throughout this paper as X). Fermiophobic bosons dominantly couple to standard model 
gauge-bosons. We consider possible strategies for distinguishing a spin-0 fermiophobic res- 
onance from a spin-1 resonance (we consider both even and odd parity). We will find that 
there are two effective complementary strategies: a search for diphoton events, and a com- 
parison of events with no spectator jets to those with one spectator jet. 

Fermiophobic bosons have been studied in a variety of contexts [THT3], particularly in 
the context of a fermiophobic Higgs boson [2]. Fermiophobic bosons can arise naturally 
in models where a hidden sector couples dominantly to exotic heavy fermions or scalars. If 
these matter fields are heavy enough, then they cannot be directly produced at the LHC. 
But if these fields are charged under standard model gauge groups, then loops of heavy 
particles mediate an effective coupling to standard model gauge bosons. One-loop diagrams 
can thus mediate mixing between standard model gauge bosons and fermiophobic bosons, 
and search strategies for this effect have been well-studied [IH3J EHD] . An alternative strategy 
is to search for triple-boson couplings, where X couples to two standard model gauge bosons. 
This strategy can be especially useful in cases where X does not mix with standard model 
gauge bosons; examples include the case where X is a scalar, or when the heavy fields 
mediating the interaction are not charged under U(l)y. 

One would expect that the most effective way to search for fermiophobic bosons with triple 
boson couplings would be through the process gg — > X — > VV, where V is an electroweak 
gauge boson. This process benefits from large production rate associated with a gluon initial 
state, and clean photon or lepton signals which are possible with an electroweak boson final 
state. This channel was studied in [12] for the case where X is a pseudovector. A variety of 
experimental searches for fermiophobic bosons have been conducted |15j . 

It is important to determine all quantum numbers of any exotic resonance discovered at 
the LHC. The measurement of the spin of a resonance has been studied by many authors, 
particularly in the context of resonances similar to the Higgs [16J. A typical method of 
spin determination is to look at the angular distribution of the decay products in the rest 
frame of the decaying particle. We will consider two complementary approaches to spin 
determination which do not require the use of angular information, but instead rely on 
effects associated with Landau- Yang theorem [T7J [18] selection rules. 

According to the Landau- Yang theorem, a massive vector boson cannot decay to two 
identical massless vector bosons. Thus, if the decay X — > 77 is observed, it is clear that 
X is not spin-1. This is the reason why it is clear that the recent exciting discovery at 
the LHC of a boson with a mass of ~ 125 GeV [19] cannot be a spin-1 resonance. The 
Landau- Yang theorem also implies that a massive vector cannot be produced from an initial 
state of two on-shell gluons. This implies that the production of a spin-1 fermiophobic boson 
from an Xgg vertex must be accompanied by the emission of a spectator jet. We will find 
that the these two strategies can be used together to distinguish between spin-0 and spin-1 
fermiophobic bosons for a wide range of models. The efficacy of each strategy depends 
primarily on the relative strength of the coupling of X to SU(2) and to U(l)y gauge bosons. 

In section II, we describe the coupling of fermiophobic bosons to standard model gauge 
bosons in terms of higher-dimension effective operators. In section III we describe the pri- 
mary X production and decay processes. In section IV we describe our method of simulating 
signal and background events, and in section V we describe the sensitivity of the LHC, and 
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its ability to distinguish the spin of fermiophobic resonances. We conclude in section VI 
with a discussion of our results. 



II. EFFECTIVE THEORY OF THE FERMIOPHOBIC BOSON 

The theory we consider consists of a new massive boson X which has effective couplings 
to standard model gauge bosons. We will assume X is not charged under standard model 
gauge groups; the coupling can thus be written as an effective operator of dimension 5 or 6 
(depending on the spin of X) which couples X to standard model field strengths. We now 
list the possible (lowest dimensional) effective operators for the four spin/parity assignments 
of this boson. 

For either a scalar or pseudoscalar, there is one lowest dimension (dimension 5) effective 
operator. 

O s = ^XTr[F, u F^}, (1) 
O ps ^^e a0jS XTr[F^F^\. (2) 

For a (pseudo)vector X, the lowest dimensional operator is of dimension 6. For the vector, 
we find four possible operators: 

Ol = ^Tr[F, a d^} (3) 

2 v = ^Tr[d^ a F^} (4) 

Ol = j- 2 d X^[F m F«e\ (5) 

Ot=j- 2 X^Tr[F a ^F a % (6) 

where we have written each operator only to quadratic order in standard model gauge field 
strengths. For a fully gauge-invariant operator, the partial derivative should be replaced by 
a covariant derivative. The first three operators are related by an integration by parts. 

X^TriF^dpF^} = dp{X»Tr[F m F a P\) - {X» Tr^F^F^} + d p X» Tr[F MQ F^]). 

Dropping the surface term, we find 

0\~ -{(% + <%). 

We will in what follows assume X to be on-shell. In this case, we find that the operator 0\ 
does not contribute. The vertex function for this operator is 

T^(k x , h, k 2 ) = + - g pv ki ■ fca) 

= ^(tfA£-<r*i-k). (7) 

Contracting this vertex function with the polarization vector for the X thus gives zero. The 
vertex function for Of, is identical, and so also does not contribute. Thus the only operators 
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which contribute to the X coupling are 0\ and 0%, which are equivalent by an integration 
by parts. We thus drop superscript labels and write 

O v = ^X»Tr[F m dpF a P}. (8) 

Finally, for the pseudovector X, the only non- vanishing operator is |T2j 

O pv = ^X^ vpa Tr[F^F^]. (9) 

For non-Abelian fields F£g, covariant derivatives should be used to act on the field strengths. 
The covariant derivatives will give rise to higher point vertices, allowing decays to 3 or more 
standard model bosons]]] For coupling to gluons, the higher point interactions have a sizable 
contribution to the production channel. 

We will assume a higher-dimensional effective coupling to all gauge groups of the standard 
model, though for generality we will allow the couplings to be different. We will thus 
characterize the couplings by an energy scale A, and encode the relative strength of the 
coupling to SU(2)l and U(l)y in relation to SU(3)qcd by C\ and Ci- The effective operators 
are thus given as 

o s = 1~xg; u g^ + ^-xw; u wr + y s xb ^ b " u ^ 

O ps = ±XG%,Q? + ^XW; u Wt u + ^-XB^ (11) 

l*-ps ^ps ^ps 



1 c c 

f.'ii , in n / itm.1 J. i-nn-i n TTrUX0 _|_ J^t 1 R D „ R°^ 

A 2 

V 



a = — X"G a m D G^ + -±X»W; a D p W ia ? + -^X^B^DpB^ (12) 



"-v V 



O pv = —X»Gl a D p G aa ? + -f X^W^DpW^ + -fx^B^DpB^ (13) 

pv pv pv 

with F^ u = \e tlvaT F ar . Since field strengths are parity even objects and e MJ/(JT is parity odd, 
we see that the X in Eq. (11) and Eq. (13) is a pseudoscalar and pseudovector, respectively. 



III. PRODUCTION AND DECAY CHANNELS 



Production at the LHC will be primarily through the channels gg, gq, gq, and qq. Vec- 
tor boson fusion is also a viable production channel for the X [10J, but will only become 
important when coupling to SU(3) is small or zero. Thus we will not consider VBF as a 
production channel in this analysis, though we note that inclusion of this production chan- 
nel will only add to the 2-jet cross section. We see (Fig. [T| that all but the gg partonic 
production channels require the presence of a spectator jet. If X is spin-1, the dominant 
gg X production channel is forbidden for on-shell X; thus spin-1 X is always produced 
with at least one accompanying jet. Since this dominant production channel is available 
to the spin-0 X, we should in principle be able to determine spin by simply comparing the 
zero-jet and one-jet cross sections for any observed fermiophobic resonance. We will return 

1 Extra terms arising from the covariant derivative will not affect two-body decay widths. They do, however, 
affect two-body branching fractions through their contribution to three-body widths, and so we use the 
full covariant derivatives in all of our couplings for event generation. 
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FIG. 1. Partonic production channels for the X 



TABLE I. Tree-level production cross sections (in fb) for X + 0-jet and X + 1-jet, as calculated 

by MadGraph5 for d = 10, C 2 = 0, and either A = 50 TeV (spin-0) or A = V^O TeV (spin-1). 

Cross sections for parity even and parity odd are identical for a spin-0 X. Jets are required to have 
p T > 50 GeV. 

m x I 500 GeV 1000 GeV 1500 GeV 2000 GeV 
= 8 TeV 

a(pp -> X {p)s ) 4.70 x 10 2 3.38 x 10 1 4.31 6.82 x 10" 1 

a{pp^X {p)v ) 

a(pp -> X {p)s + j) 2.71 x 10 2 2.69 x 10 1 4.03 7.04 x lO" 1 

<r(pp-+X v +j) 7.26 x 10 1 8.17 1.36 2.51 x lO" 1 

a(pp^X pv +j) 2.51 x 10 1 8.17 1.36 2.51 x 10" 1 
y/s = U TeV 

a(pp -> X {p)s ) 2.25 x 10 3 2.76 x 10 2 5.91 x 10 1 1.62 x 10 1 

<r(pp -> X (p)v ) 

a(pp -»• X {p)s + j) 1.53 x 10 3 2.53 x 10 2 6.40 x 10 1 1.94 x 10 1 

a(pp ->• X v + j) 7.48 x 10 2 1.32 x 10 2 3.72 x 10 1 1.24 x 10 1 

a(pp -»• X pv + j) 1.24 x 10 3 1.32 x 10 2 3.72 x 10 1 1.24 x 10 1 



to this possibility in a later section. Interestingly, it turns out that the matrix element for 
the gg —> Xg hard process also turns out to be zero if X is spin-1, so production occurs 
only through quark-gluon and quark-antiquark scattering. 

We list the production cross sections for the different spin/parity assignments in Table |l| 
as calculated in MadGraph5. Cross sections for production with zero and one jet are listed 
separately. We see, as expected, that the cross section for spin-1 X vanishes when there is 
no accompanying jet. We also see that production for spin-1 X is parity dependent at low 
masses, while for spin-0 X there is no noticeable dependence. 



A. Decay Channels 

Since the X is in general coupled to all gauge groups of the standard model, there is some 
non-zero branching fraction to both gluons and electroweak gauge bosons. If X decays via 
its coupling to gluons, the final state signal would lie in the multi-jet events; two or more 
jets for a spin-0 X and 4 or more jets for a spin-1 X (the matrix elements for X — > gg, ggg 
vanish). As this is a difficult experimental analysis, we focus instead on electroweak decays, 
which lead to much cleaner signals at the LHC. Furthermore, since we want to study an on- 
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shell resonance, it is necessary that we can fully reconstruct the X from its decay products, 
and so we do not look in channels that contain neutrinos in the final state. We will look in 
the channels X — > ZZ — > 41 (the so-called golden channel), X — > Zj — > 21 + 7, and X — > 77 
(the diphoton channel is forbidden for if X is spin-1). 
Scalar decay widths are given by 



F S (WW) _C\m\ / n Am 2 w (^ Am 2 w + g m^ . (M) 



mx 27r A 2 y m 2 ^ \ m 2 ^ m x 



r s (zz) 


_ (C 1 cos 2 ^ + C 2 sin 2 ^) 2 


m^ 


m x 


47T 


A 2 


r s (^7) 


(d - C 2 f cos 2 ^ sin 2 ^ 


m\ 


m x 




A 2 


r s (77) 


_ (dsi^^ + CsCOS 2 ^) 2 


m| 


m x 


47T 


A 2 



4mi /„ 4mi „ mi . 

+ 6^ (15) 



m x V m x m x 



r s (gg) _ 2m 2 x 
m x 7T A 2 



(16) 
(17) 
(18) 



Note that the functional structures are identical to those of the standard model Higgs decay. 
Pseudoscalar decay widths are given by 

T ps {WW) _ C 2 m\ f i _ 4<A 3/2 (1Q) 



m x 27r A 2 V m\ 



x 



Y ps {ZZ) 


_ (Ci 


cos 2 9 W + C 2 sin 2 6 W ) 2 


m| 


m x 




47T 


A 2 


T ps (^7) 


_ (Ci 


- C 2 ) 2 cos 2 9 W sin 2 # w 


m\ 


m x 






A 2 


r ps (77) 


_ (a 


sin 2 9 W + C 2 cos 2 9 w ) 2 


m| 


m x 




47T 


A 2 



4m| \ 3/2 

m 2 . ; 

m 



vr A 2 



) (20) 

(21) 
(22) 
(23) 



It is interesting to note that the scalar and pseudoscalar partial widths become equal as 
m z,w/ m x —> 0. This suggests that all information concerning the parity of the X is encoded 
in the longitudinal modes of the W and Z. 
Vector decay widths are given by, 



T V (WW) = Clm%m^ ( _ 4m|A 3/2 
m x 487T A 4 m 2 x \ m\ J 

V V {ZZ) _ (d cos 2 6 W + C 2 sin 2 6 w f m x m| / 4m? x 3/2 



m x 96n A 4 m 2 x \ m 

T v (Zj) (d -C 2 ) 2 cos 2 6 W sin 2 6 w m x m 2 z 



(24) 
(25) 



m x 96n A 4 m\ \ m\ J \ m 



m| \ / . 



1 + ^1-^ (26) 
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TABLE II. Partial widths (in MeV), as calculated by MadGraph5, for m x = 500 GeV, C\ = 10, 
C 2 = 0, and either A = 50 TeV (spin-0) or A = ^/TO TeV (spin-1). 



V X V 2 


WW 


ZZ 


Z 1 


77 


99 




3-body 


total 


r s 


680.2 


190.1 


128.8 


21.73 


31i 


\3 


89.63 


1142 


r 


677.2 


188.6 


128.8 


21.73 


31.? 


13 


90.00 


1138 




4.495 


1.634 


0.5762 










0.3815 


7.086 


- 1 - pv 


4.037 


1.416 


0.5762 










0.3516 


6.381 



Lastly, pseudovector decay widths are given by, 

T P v(WW) = Clm^m^ f _ 4mf^ 5/2 
mx 487r A 4 m\ \ m\ 

T pv (ZZ) _ (d cos 2 9 W + C 2 sin 2 6 w f m x m\ ( 4m| x 5/2 



mx 967r A 4 m\ \ m\ 

Yp^) _ (d - C 2 ) 2 cos 2 ^ sin 2 6 n 



m x m z i _ m ^ \ / . m 



mx 967r A 4 m x \ m 2 ^ / \ m\ 



(27) 
(28) 
(29) 



Again, we see that the parity-even and parity-odd widths become equal as mz,w/ m x — > 0. 

In Table [TlJ we list the widths (in MeV) to each channel for mx = 500 GeV, C\ = 10, 
C*2 = 0, and either A = 50 TeV (spin-0) or A = -\/T0 TeV (spin-1). We see that for all 
cases the narrow width approximation is valid. To see that the narrow width approximation 
remains valid throughout our parameter space, we plot the contours of r M /mi in A-m x 
space for both spin-0 (Fig. [2j and spin-1 (Fig. [3J. We neglect any contributions from three- 
body decays. We see that for spin-0, T tot /mx < 10% for a majority of the parameter space, 
though at large values of mx/A the narrow width approximation may break down. In the 
case of a spin-1 X, the narrow width approximation holds over a much wider region of 
parameter space. 




FIG. 2. Contours of T to t/mx for spin-0 X, for coupling only to gluons [C\ = C 2 = 0) (left) and 
for coupling to SU(2) with C\ = 10, C 2 = 0. For a vast majority of the parameter space width is 
less than 10% of the mass, and the narrow width approximation is valid. 
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FIG. 3. Contours of Ttot/mx f° r spin-1 X. Also shown is the line A = mx, at which point the 
effective field theory description is naively expected to break down. For coupling to SU(2) with 
C\ = 10 we find T to t/fnx *C 1 over a vast majority of the parameter space, and the narrow width 
approximation holds. 



We can now calculate the branching fractions, again ignoring contributions from three- 
body decays, from Eqs. (14)-(29). The branching fractions for spin-0 X depend on two 
parameters, C\ and C-i. We illustrate this dependence in Fig. [4] by plotting branching 
fractions for a certain choice of C± and as a function of CijC\. 




For the spin-1 X, the branching fraction to gg (and 77) vanishes, and so the overall scale 
of the electroweak couplings in relation to the gluon coupling drops out. Thus the remaining 
branching fractions only depend on C2/C1, and are constant with respect to overall scale A. 
These are shown in Fig. [5] 

Branching ratios in both cases are also dependent on mx, but for mx ~ 500 GeV the 
dependence becomes negligible. Similarly, differences between parity even and parity odd 
become negligible. 
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FIG. 5. X decay branching fractions for when X is spin-1 as a function of C2/C1, for mx 
500 GeV. 



TABLE III. Three-body branching fractions, as calculated in MadGraph5, for C\ = 10, C2 = 0, 
and either A = 50 TeV (spin-0) or A = ^/W TeV (spin-1). 



m x = 


500 GeV 


1000 GeV 


1500 GeV 


2000 GeV 


10 TeV 


X s 


0.078 


0.154 


0.207 


0.24 


0.436 


Xp s 


0.079 


0.154 


0.206 


0.245 


0.447 


x v 


0.054 


0.129 


0.199 


0.266 


0.818 


Xp v 


0.055 


0.129 


0.199 


0.267 


0.816 



B. Three-body decays 

In the previous section, we ignored three-body decay widths due to the usual assumption 
that the reduction in phase space and extra coupling constant factors render these negligible. 



There are cases, however, where this assumption cannot be made. We show in Table III 
the three-body branching ratios for all four spin/parity assignments at five different masses. 
We see that when X is spin-0, the three-body branching fractions are small, but in general 
they are not negligible. More importantly, we see that if X is spin-1, the three-body decays 
strongly dominate for mx 3> mw,z- This behavior arises again from the Landau- Yang 
theorem. For mx 3> m\y,z, the electroweak gauge bosons can be considered approximately 
massless, and the two-body decay becomes highly suppressed. We see this clearly in the 
vector and pseudovector three-body branching fractions for mx = 10 TeV. 

It is thus important to note that searches for two-body resonances will be ineffective 
for very massive fermiophobic spin-1 particles. Furthermore, any search for a three-body 
resonance will be complicated by the fact that, for electroweak decay products, full recon- 
struction using only leptons and photons is impossible. The only three-body electroweak 
decay modes for the X are W + W~ Z and W + W~ r y; purely leptonic decays will always pro- 
duce missing momentum via neutrinos. Thus for full reconstruction we must rely on the 
hadronic channels, (either X — > gqq or hadronic decays of electroweak gauge bosons). Since 
a vector X must be produced with an associated jet, this signal will be contained within 
events containing 4 or more jets. Alternatively, a search can be performed in channels in 
which intermediate electroweak bosons are off-shell. 

For m x < 2000 GeV, the branching fractions to electroweak gauge bosons dominate, and 
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we may still be able to achieve discovery in the two-body decay channels. The remainder 
of this paper will focus on this possibility, and the possibility of distinguishing spin if a 
resonance is indeed observed. 



IV. COLLIDER SIMULATION 



As outlined above, we will generate event samples in which X is produced through the 
Xgg vertex and subsequently decays through electroweak couplings. We generate parton- 
level events in MadGraph5 |20j , shower and hadronize the events in Pythia [2TJ , and perform 
detector simulation using PGS4 [22] using the ATLAS detector card. To generate the model 
files for MadGraph, we have used the Mathematica package FeynRules (23], to which we 
feed the effective interaction operators, Eq. ( 10 )-( 13 ), as input. We generated events for 
collider energies of y/s = 8 TeV and y/s = 14 TeV. 

For spin-0 X, we generate both pp — > X and pp — > Xj at the parton level, in order to 
include all possible production channels. These two processes are matched at the Pythia level 
using the MLM algorithm in order to avoid over counting between matrix element generated 
and ISR generated jet-containing events. For the spin-1 cases, events are generated solely 
from the 1-jet matrix element (which in this case is tree-level). 



A. Cuts 



The following reconstruction cuts will be applied for the three channels we study. These 
cuts are applied at the detector level. 



• X -> ZZ 

We require events to contain 4 charged leptons (e , /i 1 * 1 ) with 77 < 2.5. These must 
consist of two pairs of same-flavor opposite-sign leptons which each have invariant 
masses in the range 80 GeV < rrii+i- < 100 GeV. The total invariant mass must lie 
within 10% of m X - 

• X -> Z7 

We require events to contain a photon as well as 2 same-flavor opposite-sign leptons 
(77 < 2.5) whose invariant mass lies in the range 80 GeV < rrii+i- < 100 GeV. The 
two leptons and photon then must have a total invariant mass which lies within 10% 
of m X - 

• X — > 77 

We require events to contain 2 photons (77 < 2.5) which reconstruct to invariant mass 
within 10% of mx- 



In addition to these cuts, a number of standard isolation cuts are applied at the MadGraph 
level, including requiring AR > 0.4 for any pair of particles. At this point we do not impose 
any cuts on the jet structure of the events. 
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B. Background 



We generate background events using the same MadGraph5 simulation chain. For the 
vast majority of events where a lepton pair has an invariant mass between 80 — 100 GeV, the 
leptons arise from the decay of a Z boson. Other sources of lepton pairs satisfying the cuts 
(including lepton misidentification) are insignificant and can be ignored. For the ZZ and Z7 
background channels, we use 0-jet and 1-jet matrix elements, which are matched using the 
MLM algorithm. For the 77 channel, we use only the 0-jet matrix element. Cross sections 
for these channels, after imposing the above mentioned cuts in each X-mass window, are 
given in Table IV at LHC energies of 8 TeV and 14 TeV (see also [2U 125]). 



TABLE IV. Inclusive cross sections (in fb) for standard model background in the ZZ, Z7, and 77 
channels, after reconstruction cuts in each total invariant mass window. 













m x 












500 GeV 


1000 GeV 


1500 GeV 


2000 GeV 






» ZZ - 


-> l+l- 


-1+1-) 


2.09 x 10- 1 


7.35 x 10- 3 


1.87 x 10" 4 


2.16 x 10~ 5 


8 TeV 


<?bg{pP ~ 


* Z7 - 


> l+r 


7) 


3.93 


2.98 x 10" 1 


2.55 x 10~ 2 


2.93 x 10~ 3 




<?bg{pP ~ 


» 77) 






1.80 x 10 1 


1.42 


2.08 x 10- 1 


4.66 x 10" 2 




ctbg(pP ~ 


» ZZ - 


-> i+v 


-i+r) 


4.41 x 10" 1 


2.45 x 10" 2 


1.90 x 10~ 3 


1.61 x 10" 4 


14 TeV 


<?bg(pP ~ 


> Z7 — 


> i+r 


7) 


8.42 


4.16 x 10" 1 


1.09 x 10" 1 


2.72 x 10~ 2 




ctbg(pP ~ 


+ 77) 






3.65 x 10 1 


2.79 


8.73 x 10" 1 


3.27 x 10" 1 



V. COLLIDER REACH 

Since the narrow width approximation holds, the number of signal events in a given 
channel and at a given luminosity £ depends on the coupling constants G\ and C 2 only 
though the final state branching fractions. If we denote by cr pro d the cross-section for the 
production process pp — > X + (0 and/or 1 — jet), and by e the fraction of such events which 
pass the analysis cuts in any channel, then the number of signal events is given by 

N sig = a pmd x e x BR(X -> VV) x C 

oc AT n x BR(X -»■ VV) x C (30) 

where n = 2 for spin-0, and n = 4 for spin-1. For a discovery, we will require a 5a (Gaussian 
equivalent) excess of signal events over background events satisfying the same cuts, and we 



require at least 5 signal events. With Eq. (30) in mind, we define the LHC reach for any 



integrated luminosity in terms of the quantity 

R vv = A/[BR(X -> VV)] l/n . (31) 

We find the reach for 10 fb -1 and 30 fb -1 at center of mass energy of y/s = 8 TeV, 
and 100 fb -1 at y/s = 14 TeV. For the scalar and vector cases, the LHC reach is plotted 
against mx in Fig. [6] We plot the inclusive X — > ZZ channel (both and 1-jet events) 
and the inclusive X — > Z^ channels separately. The reach plots for the pseudoscalar and 
pseudovector case are essentially identical to the scalar and vector cases, respectively. We 
see that in both channels, we can probe higher scales in the spin-0 case, which is expected 
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FIG. 6. Collider reach for both ZZ (top) and Z7 (bottom) channels. Spin-0 X is shown at left, 
spin-1 X at right. The reach plots for the parity odd bosons are essentially identical to their parity 
even counterparts. 



from the additional A 2 suppression coming from the dimension 6 operators coupling a spin-1 
fermiophobic boson to the standard model. We also note that the Z7 channel is generally 
a better probe despite the higher backgrounds. This is mainly due to the small Z — > l + l~ 
branching fraction and the reconstruction cuts present in the ZZ channel. 

In addition, we show in Fig. [7] the reach for the 77 channel, present only for spin-0 X. 
Current LHC data can probe well into the hundreds-of-TeV range. 
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FIG. 7. Collider reach for both 77 channel, which is present only for spin-0 X. 



Note that we have not accounted for NLO QCD corrections to the X production cross 
section. The associated ^-factors are usually greater than one, indicating that a correct 
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treatment of NLO effects would likely increase the LHC reach. 



A. Spin Determination 



It is clear from the reach plots of the previous section that the simplest way to distinguish 
a spin-0 resonance from spin-1 is to look to the 77 channel. If there is an observation in the 
77 channel, with or without a corresponding observation in ZZ or Z7, the Landau- Yang 
theorem implies that the decaying particle is not spin-1. 

There are cases, however, where a spin-0 X will not decay to two photons. We see 



from Fig. [4] and Eq. (17) that the branching ratio BR(X — > 77) drops to zero for 
C2/C1 = — ta,n 2 8 W (an explicit model where this occurs is exhibited in the Appendix). 
The electroweak decay channels for the spin-0 and spin-1 X are in this case identical. We 
must in these cases resort to an alternative method of determining spin. 

Another way of distinguishing between spin-0 and spin-1 resonances is from the presence 
of extra jets; the 0-jet process pp — > X — > ZZ, Z7 is possible if X is spin-0, but not if it 
is spin-1. As we saw, we can define the LHC reach for X (spin-0) in terms of the quantity 
R vv = A/BR(X VV) 1 ' 2 . We thus see that a search in the ZZ or Z7 channel with zero 
extra jets will be more promising than a search in the 77 channel if 



R 



zz,z~ t 



(0-jet) 



R 



77 



BR{X ZZ, Zi) 
BR(X -> 77) 



> 1. 



(32) 



The ratios BR(X -> ZZ)/BR(X 77) and BR(X Zj)/BR(X 77) depend only 
on the quantities m x and C 2 /Ci. Having determined Rzzi^x) and Rz^imx) for the 0-jet 
sample from the generated signal and background events, one can determine the range of 
C2/C1 (for any mx) over which the ZZ or Z7 0-jet search provides better prospects for 
discovery than the 77 channel. These ranges are plotted in Fig. [8j note that these ranges 
would not change significantly after the LHC energy and luminosity upgrade. 
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FIG. 8. Range of C2/C1 for which a ZZ-channel search (left) or a Z7-channel (right) search is 
more promising than a diphoton search. 

Thus far, our analysis has focused on strategies for determining if X is not spin-1. One 
might ask the converse question: does the appearance of an excess in ZZ and/or Zj events 
with one extra jet, unaccompanied by an excess in 0-jet events or 77 events, necessarily 
imply the X resonance is spin-1? Essentially, this amounts to the question of whether or 
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not a spin-0 resonance can produce an excess in 1-jet events, without also producing a 
statistically significant excess in 0-jet events or 77 events. 

We can determine this by comparing the LHC reach for 0-jet events to the reach for 1-jet 
events in the ZZ or Z7 events, assuming X is spin-0. If the reaches are comparable, then a 
spin-0 X coupling to the standard model in such a way as to produce a 5a excess in ZZ or 
Z7 events with 1 extra spectator jet would, for similar luminosity, also produce an excess 
in events with no spectator jets. An observed 5a excess in 1-jet events without some excess 
in 0-jet events would imply that X is not spin-0. 

In Fig. |9]we plot Ryy (0-jet) / Ryy (I- jet) as a function of m x for the ZZ and Zj channels 
(assuming X is spin-0). From this plot, we see that for much of the range, an observation 
of a resonance with one spectator jet, but an absence of a excess with zero spectator jets, 
is sufficient to demonstrate that the resonance is not spin-0. Furthermore, we see that 
this conclusion is largely independent of the collider energy, luminosity, and decay channel 
studied. We also show (Fig. 10) Ryy (0-jet) / Ryy( 1-jet) for the case of a spin-1 X (in this 
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FIG. 9. Ratio of reaches Ryy (0-jet) /Ryy( 1-jet) for X ->■ ZZ (left) and X -> Z7 (right) in the 
case of spin-0 X. We see that if the resonance is spin-0, we expect a similar number of 0-jet and 
1-jet events. 



case, the number of signal events is proportional to R~ A ). As expected, we see that an 
observation in the 0-jet channel without a corresponding observation in the 1-jet channel is 
sufficient to demonstrate that the resonance is not spin-1. 



VI. CONCLUSIONS 

The existence of exotic fermiophobic bosons is well-motivated theoretically. These types 
of bosons could reveal themselves as resonances which decay to two or more SM bosons at 
the LHC. We have established the efficacy of the LHC in probing models of this type, and 
demonstrated the viability of two methods for determining spin. The existence of a diphoton 
decay channel establishes that the resonance is not spin-1, as is well-known. When a boson 
does not decay to two photons, spin can be determined through an inspection of the jet 
structure of the signal events. Over much of the parameter space, a discovery and a spin 
determination can be achieved in as few as 5 events, thus requiring far fewer statistics than 
alternative methods of spin determination (e.g. analyses of angular distributions of decay 
products). 
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FIG. 10. Ratio of reaches Eyy(0-jet)/ifyy(l-jet) for X ^ ZZ (left) and X -> Z7 (right) in the 
case of spin-1 X. We see that if the resonance is spin-1, we expect to many 1-jet events and very 
few 0-jet events. 



The LHC has already discovered a new boson, and the existence of diphoton decays shows 
definitively that this new particle is not spin-1. We have shown that the LHC has very good 
prospects for discovering fermiophobic bosons in the near future, even if they are relatively 
heavy, or are only coupled to the standard model by higher dimension effective operators 
which are heavily suppressed. For example, a 100 fb _1 run of the LHC at 14 TeV could 
find 5<T evidence for a scalar with mx ~ 2 TeV, coupled to the standard model by effective 
operators suppressed by a mass scale ~ 300 TeV. The reason for this large reach is that a 
fermiophobic boson can be produced from gluon couplings, but observed through electroweak 
decays, a channel which is ideal for detection at the LHC. In particular, the diphoton channel 
is clearly the most promising, as the signal is very clean and the background is very small. 

But one typically expects an electrically neutral boson to have a very small branching 
fraction for decay to 77. This is the case for the 125 GeV boson discovered at the LHC. 
Fermiophobic scalars are a major contrary example; the branching fraction for diphoton 
decay can easily be 0(1). One can see why such large X — > 77 branching fractions are 
allowed by considering a high-energy theory where the fermiophobic boson X couples to 
standard model gauge bosons only through loops of heavy fermions and scalars which are 
charged under standard model gauge groups. In this scenario, since all decays arise from 
one-loop diagrams, the X — > 77 branching fraction can be comparable to that of other 
channels. This behavior is markedly different from other scenarios in which some decays to 
standard model particles occur at tree-level. Fermiophobic bosons are thus one of the most 
interesting prospects for future discovery at an upgraded LHC. 

Although the methods we have discussed here are applicable for spin determination, it 
is clear from our results that parity information is generally not encoded in jet number 
distributions. Alternative methods must be used to determine the parity of the couplings, 
which would give clues to the nature of the underlying physics of the new resonance. 
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Appendix A: Vanishing BR(X — > 77) 

The condition that C 2 = —C\ tan 2 9 W , which is necessary for a vanishing branching ratio to 
two photons, can arise naturally as follows. Let us assume that the X couples to electroweak 
gauge bosons through triangle diagrams involving extra heavy fermions. In this model there 
are two extra vector-like fermion multiplets. One is an SU(2) doublet, 



, with Y = -1, 



which couples to X through a Yukawa interaction with coupling constant Ai. The other is 
an SU(2) singlet, c, with Y = —2 and which couples to X through a Yukawa interaction 
with coupling constant A2. The electric charge (Q = T 3 + Y/2) of b and c is — e while the 
charge of the a is zero. 

Let us first consider the coupling to two photons. Only b and c will contribute to the 
loop, so the amplitude is proportional to 




EM 



EM 



oc Ai(-e) 2 + A 2 (-e) 



We see that for the X77 coupling to vanish, we must have Ai = — A 2 . Now, we look to the 
couplings to the SU(2) and U(l)y gauge bosons. For SU(2), only the doublet runs in the 
loop, and we have 

■ 517(2) 

oc Ai# 2 Tr(titj) 
SU(2) cx|A l5 2 

We can identify this factor as our coefficient C\. For U(l)y, all three particles flow in the 
loop, and we find 





ocA^ 2 (-|) 2 + A^ 2 (-I) 2 + W 2 (-l) 2 
= \X ig ' 2 + W 2 



1 \ 
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where in the last line we replaced A2 = — Ai. This coefficient can be identified with C 2 . 
Since g' = gtan9 w , we see that C 2 = —C\ tan 2 6 W as needed. 

The above argument is strictly only valid if the particles coupling to the X — a, b, and c 
- have identical masses. If this condition is relaxed, then the diphoton coupling will vanish 
only for Ai ~ — A 2 . 
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